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Majorana Fermion Rides on a Magnetic Domain Wall 
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We propose using a mobile magnetic domain wall as a host of zero-energy Majorana fermions in 
a spin-orbit coupled nanowire sandwiched by two ferromagnetic strips and deposited on an s-wave 
superconductor. The ability to control domain walls by thermal means allows to braid Majorana 
fermions nonintrusively, which obey non-Abelian statistics. The analytical solutions of Majorana 
fermions inside domain walls are obtained in the strong spin-orbit regime. 
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Introduction. —Majorana fermions (MFs) have at¬ 
tracted much attention due to their non-Abelian ex¬ 
change statistics, which makes them potentially useful 
for topological quantum computation In a film of 
a Rashba spin-orbit coupled (SOC) semiconductor sub¬ 
jected to a Zeeman field and proximity-induced s-wave 
superconductivity, a vortex in the superconducting order 
parameter supports one MF at its core ma, as it does in 
a, p + ip superconductor [3] . A one-dimensional wire with 
the same ingredients also supports MFs at the bound¬ 
ary between topological and nontopological regions. Such 
wires can form a mesh, where the MFs can be braided 
by slowly adjusting gate voltages, as proposed by Alicea 
et al. la. 

Spintronics aims at the active control and manipula¬ 
tion of spin degrees of freedom in condensed-matter sys¬ 
tems [a. Topologically stable magnetic textures, e.g., a 
domain wall (DW) in an easy-axis ferromagnetic wire or 
a vortex in an easy-plane ferromagnetic film, have been 
extensively studied out of fundamental interest as well as 
practical motivations exemplified by the racetrack mem¬ 
ory 0. Their dynamics can be driven by various means, 
e.g., an external magnetic field [S], an electric current 
(in metallic ferromagnets) [^ , or a temperature gradient 
mm- The topological phase transition of the SOC 
nanowire with proximity-induced superconductivity can 
be driven by adjusting the effective magnetic field na¬ 
if we can identify a localized magnetic texture that sup¬ 
ports MFs under suitable conditions, we would, therefore, 
be able to manipulate them by controlling the magnetic 
host with standard spintronic techniques. 

In this Letter, we propose to use a mobile DW in a 
proximate ferromagnet as a MF host. We show how to 
exchange MFs in the Y junction by thermally-driven mo¬ 
tion of the DWs, exhibiting non-Abelian braiding statis¬ 
tics. In conjunction with ongoing rapid developments 
in control and motion of magnetic DWs, it could be 
envisioned that networks of magnetic wires in proxim¬ 
ity to a superconducting substrate can serve as robust 
and versatile means for storing and manipulating MFs. 
Specifically, we obtain analytic solutions for MFs inside 
a smooth DW in the strong spin-orbit regime. An intu¬ 



FIG. 1. (Color online) A schematic diagram of a device sup¬ 
porting two MFs in a SOC nanowire with proximity-induced 
exchange field and superconductivity. The wire is sandwiched 
between two easy-axis ferromagnets: one with a uniform mag¬ 
netization and the other with a DW. In the right portion of 
the wire, where the magnetizations of the two ferromagnets 
are parallel, the net exchange field engenders the topological 
region, fomenting MFs at locations marked by A and B. The 
temperature-gradient induces magnon drag in the x direction, 
which drives the DW and its accompanying MF at A to the 
hotter region. 


itive idea of the existence of MFs inside the DW can be 
developed by casting the Bogoliubov-de Gennes (BdG) 
equation for MFs as two Jackiw-Rebbi Dirac equations 
with spatially-varying mass terms m- 

MF at a magnetic DW. —Consider a semiconducting 
nanowire with Rashba SOC, which is deposited on an 
s-wave superconductor and is proximity-coupled to ad¬ 
jacent ferromagnets [a [3]. See Fig. for a schematic 
design of a device. The superconducting quasiparti¬ 
cle spectrum is obtained by solving the BdG equation 
where 

dl - p + iad^u^ Ta -I- M • (T -f Ari (1) 

acts on the spinor wavefunction T' = . 

Here, m, p, and a are the electron effective mass, 
the chemical potential, and the strength of SOC. The 
proximity-induced exchange field M is perpendicular to 
the spin-orbit field (x y. The proximity-induced super¬ 
conducting order parameter A is gauge shifted to be real 
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and positive. Pauli-matrix vectors cr and t act respec¬ 
tively on the spin and the electron-hole subspaces of the 
spinor The corresponding quasiparticle creation oper¬ 
ator is y'l' = / dxY,a=^^]\Ua{x)'li)l^{x)+Va{x)'4ja{x)]. For 
the uniform exchange field, the “topological gap” at zero 
momentum is given by 

Eg = |M| - (2) 

When the gap Eg is positive, |M| > the wire 

is in the topological phase, harboring a pair of MFs at its 
ends [121 El- Otherwise, for the negative gap, the wire 
is in the normal phase, without MFs. 

A spatially-varying exchange field induces the topo¬ 
logical phase transition along the wire where |M| crosses 
y/A2 -f ^2 A DW in a ferromagnet adjacent to the 
wire is a natural object to bring about such a position- 
dependent field. We assume that the energy of the ferro¬ 
magnet is given by 

t/[m(a;)] = J dx - Ka:ml + Kyinl]/2, (3) 

where m(a;) is the unit vector in the direction of local 
magnetization. Here, A, Kj., and Ky are the positive 
coefficients characterizing the stiffness of the magnetiza¬ 
tion, the easy-axis anisotropy along the wire, and the 
hard-axis anisotropy in the spin-orbit field direction, re¬ 
spectively |8]. Two ground states have uniform mag¬ 
netization m = ±x. A DW is a topological soliton 
solution minimizing the energy (§ with the boundary 
condition m(a; —)■ ±oo) = ±x: its magnetization is 
m(a:) = tanh(a;/A)x -|- sech(a;/A)z where A = A/K^ 
is the DW width [5]. 

We sandwich the wire between two ferromagnets: one 
with a DW and the other with a uniform magnetization, 
as shown in Fig.j^ The proximity-induced exchange field 
is described by 

M(a;) = Ml [tanh(a;/A)x -I- sech(a;/A)z] -|- M 2 X, (4) 

which introduces a new length scale A to the Hamiltonian. 
The presence of the DW causes spatial variance of the 
gap: 

^gi^) = \J + ^2 + 2 M 1 M 2 tanh(a:/A) — a/A2 -|- fjp. 

_ (3 

For example, when Mi = M 2 and 2Mi > A'^ + fE, 

there must be a pair of MFs in the wire: one at the 
right end, x —>■ -l-oo, and the other one at the DW, Xq = 
Atanh“^[(A^ -I- /i^)/2Mf — 1]. The topological stability 
of the DW protects the hosted MF from disturbances of 
magnetization. 

The technique to control a DW has been developed 
over decades in spintronics [B], which is directly trans¬ 
lated into the ability to manipulate MFs. DWs are con¬ 
ventionally driven by an external magnetic field [5] or 



FIG. 2. (Color online) The exchange of two MFs bridged by 
a topological region performed by controlling the exchange 
field, (a)-(g) The positions of MFs can be exchanged by a 
series of DW displacements, (h), (i) Uniform rotation of the 
magnetization between MFs by 180° about the spin-orbit field 
direction transforms the state (g) to the initial state (i). This 
is essentially monodomain flipping. 


a spin-polarized electrical current (in an itinerant ferro¬ 
magnet) [3]. These methods, however, affect the elec¬ 
trons in the wire by altering the Hamiltonian. Instead, 
we propose to induce the motion of a DW by thermal drag 
[10] . which is not intrusive to the electrons (see Fig.[^. A 
temperature gradient generates a magnon current flowing 
to the colder region. Magnons adjust their spins toward 
the local direction of magnetization and react by exerting 
a torque on the DW. Addition of angular momentum to 
the DW causes its shift toward the hotter region. The 
resultant velocity of the DW, therefore, is proportional 
to temperature gradient. Multiple DWs on the single 
wire can be moved simultaneously, thereby providing the 
ability to control a series of MFs. 

Exchange of two MFs .—The exchange of two MFs is 
possible in the Y junction of three nanowires, each of 
which is sandwiched between two ferromagnets. Figurej^ 
shows the process to exchange MFs bridged by a topo¬ 
logical region m- This proposal uses a well-controllable 
topological soliton (DW) as a host of MFs, and thus 
does not need intricate gate fabrication and control, un¬ 
like the exchange process proposed by Alicea et al. |S]. 
The non-Abelian exchange statistics of MFs is a universal 
property, meaning that it is invariant under the contin¬ 
uous deformation of the Hamiltonian as long as the gap 
Eg{x) r emains finite, except at isolated locations of 
MFs |5fn6|. Our proposal uses the spatially varying ex¬ 
change field M(r) and the chemical potential fixed at 
zero, /r(r) = 0. These two functions, M(r) and fr(r), can 
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FIG. 3. (Color online) (a) The trajectories of the parameters 
{cp, ri) of MF A and B in the proposed exchange process shown 
in Fig. A circle shows the state of each step in the figure. 
The combined path of two trajectories forms a closed cnrve 
on the parameter space which is topologically a torns. (b) 
The trajectories of {(f>, r]) in the process proposed by Alicea 
et al. [^. Two combined paths of (a) and (b) are topologically 
equivalent. 

be continuously deformed to a uniform exchange field in 
the z direction and a spatially-varying chemical potential 
of the form studied by Alicea et al. , without changing 
the gap. This is accomplished by first locally rotating the 
exchange field about the spin-orbit field direction to be 
oriented along the z axis: M'(r) = M(r)z, while keeping 
the gap Eg(x) unchanged. The following transformation 
then connects our Hamiltonian (s = 0) to that of Ref. [S] 
(s = l): 

M (r; s) = M (r) -|- [— M(r) -|- Mg] s, 

/- ( 6 ) 

^(r; s) = Y {[-M(r) -|- Mg] s -I- A}^ - A^. 

where Mg = maxrM(r). The gap-preserving transfor¬ 
mation between two exchange processes, one performed 
by controlling the exchange field and the other by gating 
the chemical potential, ensures that the process shown 
in Fig. [^exhibits the non-Abelian exchange statistics ob¬ 
tained in Ref. [5]. 

Alternatively, we may understand the topological 
equivalence of two processes following the general ap¬ 
proach developed by Halperin et al. [TS]. The state 
of each MF, denoted by A and B, can be described 
by two parameters. One is the wire orientation w“ = 
(cos(^“, sin(/)“) (a = A, B) in the ccy-plane pointing from 
a nontopological region to a topological region, e.g., = 

ll 7 r /6 and 4>^ = 57 r /6 in Fig. [^a). The other parameter 
is the exchange-field direction b“ = cos 77“ z -|- sin 77“ w“ 
acting on a MF, e.g., 77 ^ = — 7 r /2 and 77 '® = 7 r /2 in 
Fig. [^a). The state {(jjjrj) of each MF evolves under 
the exchange process while making a trajectory on the 
parameter space which is topologically a torus. When 
the Hamiltonian goes back to the original form after the 
exchange, the combined path of trajectories of two MFs 
forms a closed curve on the torus. Figure [^a) and (b) 
depict the combined curves of the parameters of our pro¬ 
cess and Alicea et al. [^’s process, respectively. Both 


curves wind the torus once along the angle (j): they are 
topologically equivalent. Thus the two processes should 
exhibit the same exchange statistics [16]. 

Exact solutions for MFs. —In the strong SOC regime. 
Ego = mo? max(|M|, A), we can use the BdG 

Hamiltonian density linearized at zero momentum: 

= -iad^dir-i + + An, (7) 

where the chemical potential is set to zero mini. In 
the following discussions, we set a = 1. 

To obtain an intuitive idea of the existence of MFs 
bound to a DW, let us neglect the exchange field 
for the moment. The Hamiltonian density 77^“, then, 
can be classified as the BDI class with the particle-hole 
symmetry operator a 2 T 2 K, the BDI time reversal sym¬ 
metry operator aiXiK (different from the conventional 
one), and the BDI chiral symmetry operator [IB] . 
The Hamiltonian density 77*“ and the BDI time rever¬ 
sal symmetry operator can be block-diagonalized simul¬ 
taneously by employing Majorana operators instead of 
fermion operators. Define four Hermitian Majorana op- 
erators 7 -^ ^ in such a way that 

+ ^7f )/v^> = ( 7 ® + *7^)/v^- (8) 

Majorana operators are labeled by A and B for the con¬ 
venience of the following discussion. The quasiparticle 
creation operator is 7+ = / dxY,a=t,i > 

where the new 4-component spinor ^ = 

is related to the original spinor 
by the unitary transformation d' = [/'F. In this new 
basis, the 4-component Schrodinger equation 19(4' = 
([/77*“C/'^)^ is cast as two Dirac equations 

+ [—Mx{x) — A]u^ = 0, (9a) 

i'y^dgU^ + [Mg;{x) — A]u^ = 0, (9b) 

where we have taken the representation of Dirac ma¬ 
trices to be 7° = —<72 and 7* = (— for u^). 

The masses of A and B, m^{x) = —M^{x) — A and 
m^{x) = Mx{x) — A, are essentially the topological gap: 
Eg{x) = m^{x) when Mx{x) > 0 and Eg{x) = m^{x) 
when Mx{x) < 0. These Dirac equations coincide with 
that of solitons with fermion number 1/2 that Jackiw 
and Rebbi [13] studied. Each Dirac equation always 
has a unique static normalizable solution when its mass 
changes sign between two ends a; —7- ± 00 . The exchange 
field M = [Ml tanh(a;/A) + M 2 ]x yields two zero-energy 
solutions [19] 

7 -® = y dxe~^’^^~^^^sech^^^{x/X)f^{x) 


(10b) 
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teristic degree of our problem. Figure shows the pro¬ 
file of the topological gap and the amplitude squared 
of the MF solution. The relative magnitude of the 
two components can be obtained in two limiting cases: 
\u^{x)\/\u^{x)\ ~ 2MiX when MiA 1 , and ~ 1 when 
Ml A ^ 1. For an abrupt wall A —0, the solution 0 
converges to the solution ( 10 a) which was obtained in the 
absence of M,. 


FIG. 4. (Color online) (a) The topological gap Eg{x) for 
the exchange field M = [2-|-2 tanh(a:/A)]x-I-sech(a;/A)z with 
domain wall widths A = 0.2, 1, and 5. (b) The ampli tud e 
squared |4'(a;)p = -|- |u:f( 2 ;)|^ of MFs in Eq. (Ill, 

which are located most likely where the gap vanishes. 


up to the normalization factor, and 7 ^ are normaliz¬ 
able MF solutions when M 2 < Mi — IS. and M 2 < Mid-A, 
respectively. 

When the uniform held M 2 is smaller than Mi — 
A, both masses m^{x) and m^{x) cross zero at 
Atanh“^[—(M 2 ± A)/Mij (- 1 - for A) and thus two MFs 
exist [5D]. Perturbations breaking the BDI time re¬ 
versal symmetry split the degeneracy. First, the ki¬ 
netic energy {p?d^/2w^ omitted from the original BdG 
Hamiltonian density (iH hybridizes MFs with the energy 
< max(M^, A^)/2TOa^, which is hnite, yet still much 
smaller than the bulk gap ~ Mi, A. Second, note that 
both MF operators are made of f-spin fermion operators 
|21j . The physical fermion / = ( 7 '^-|-i 7^)/2 constructed 
from them, thus, couples to the exchange held in the 2 ; 
direction. For M^{x) = M^sech{x/X), the coupling en¬ 
ergy of the fermion / to the exchange held, i.e., the hy¬ 
bridization energy of MFs, is proportional to M^A. An 
inhnitely abrupt wall A —>■ 0 reproduces the case of the 
unidirectional exchange held. 

When the uniform held exceeds Mi — A, while still 
smaller than Mi -|- A, the left part of the wire x 
— A is driven into the nontopological phase, destroy¬ 
ing MF A. The surviving MF B is protected from per¬ 
turbations, e.g., hnite M^, as long as the bulk gap 
remains hnite. In particular, for the exchange held 
M(x) (|^ with Ml = M 2 > A/2, one MF is located 
at xq = Atanh“^(A^/2M/ — 1). The anticommutativity 
of the particle-hole transformation operator a 2 T 2 K and 
the Hamiltonian density 0 allows us to obtain a 
nondegenerate MF solution given by 


„Aa: 


(lla) 


+ (,Vl + e-2AA 

uf{x) = -{v + \I2)-^^=P-YJ^ f , ^ ^ 

M ^ ^ ^ + e2^/^ VvTTi=WAy 

(lib) 


up to the common normalization factor, where is 

the associated Legendre function of the hrst kind with 
degree a and order /? and v = 2 MiA — 1/2 is the charac¬ 


Discussion .—We have proposed to bind MFs to mag¬ 
netic DWs in the heterostructure of a SOC nanowire, an 
s-wave superconductor, and ferromagnet wires. We have 
also studied MFs analytically in the strong spin-orbit 
regime. A topologically stable DW provides a robust 
matrix for strong MFs. Typical strong SOC semicon¬ 
ducting nanowires (e.g., InAs), magnetic insulators (e.g., 
EuO), and s-wave superconductors (e.g., Nb) provide the 
parameters a ^ 5meV nm, |M| ^ ImeV, A ^ 0.5meV 
m, which must be sufficient to bind MFs to DWs. Static 
MFs at DWs can be experimentally observed by mea¬ 
suring the zero-bias peak | 22 ] . free from spurious “end” 
effects such as localized states near the wire’s end. 

Diverse means to control DW motion allow us to ma¬ 
nipulate MFs. In particular, we have proposed a pro¬ 
cess braiding two MFs in the Y junction performed 
by thermally-driven DW motion, which exhibits non- 
Abelian exchange statistics. Thermally-driven motion of 
a DW has been observed in yttrium iron garnet films 
m- the DW moves at the velocity v ~ 100 /im/s for a 
temperature gradient VT ^ 2^eV/^m. The resultant 
temperature drop over the DW width A ^ 60 nm is much 
smaller than the induced topological gap ~ 200 /reV [12] . 

A domain wall in an easy-axis magnetic wire is a 
one-dimensional topological soliton. It would be worth 
pursuing higher dimensional generalizations of our pro¬ 
posal. For example, a thin film of Feo. 5 Coo. 5 Si subject 
to an external held supports a lattice of skyrmions which 
are two-dimensional topological solitons |23|. For a sin¬ 
gle skyrmion, the local direction of the magnetization 
changes smoothly from -|-z at its center to —z at its out¬ 
side (or vice versa). The magnetization on its radial line 
dehnes a domain wall located at its edge where the mag¬ 
netization is in-plane. Skyrmions, therefore, may host 
MFs on their edges under suitable conditions, which can 
be braided by thermally-induced motions of skyrmions 

m- 
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